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We study the effect of coupling magnetic impurities to the honeycomb lattice spin— 1/2 Kitaev 
model in its spin liquid phase. We show that a spin-S* impurity coupled to the Kitaev model 
is associated with an unusual Kondo effect with an intermediate coupling unstable fixed point 
K c ~ J/S separating topologically distinct sectors of the Kitaev model. We also show that the 
massless spinons in the spin liquid mediate an interaction of the form Sf a Sjp/Rij between distant 
impurities unlike the usual dipolar RKKY interaction Si a Sj a / R% noted in various 2D impurity 
problems with a pseudogapped density of states of the spin bath. Furthermore, this long-range 
interaction is possible only if the impurities (a) couple to more than one neighboring spin on the 
host lattice and (b) the impurity spin S 7^ 1/2. 



Impurity effects are an essential part of our under- 
standing of strongly correlated electron systems , both as 
a probe for the underlying electronic state as well as due 
to the numerous nontrivial effects they have on the prop- 
erties of the system 0, 0] ■ Recently many studies have 
been made of impurity effects as a probe for the putative 
quantum spin liquid state in underdoped cuprate super- 
conductors 3 and geometrically frustrated magnets 
The 5 = 1/2 honeycomb lattice Kitaev model [f| pro- 
vides a very appealing playground in this context - it has 
a gapless spin liquid phase and short range spin correla- 
tions [6] making it different from many other extensively 
studied spin liquids @, 0, 0); and crucially, the model 
is integrable via several schemes of spin-fractionalization 
into fermions 

SI- 

The Kitaev model has been studied 
in various contexts ranging from thepossibility of quan- 
tum computation with the anyons [5|, [9( that the model 
predicts, understan ding dynamics of quantum quenches 
in a critical region [10( to fractional charge excitations 
in topological insulators [Til] - However no study of mag- 
netic impurity effects in the Kitaev spin liquid has yet 
been made, which constitutes the subject of this paper. 

We study the behavior of spin-5 impurities in the gap- 
less spin liquid regime of the Kitaev model on the hon- 
eycomb lattice. The impurity coupling K scales away 
from an unstable fixed point K c ~ J/S irrespective of 
the sign of impurity coupling, similar to impurity prob- 
lems in pseudogapped bosonic spin liquids [l2j. The Ki- 
taev magnetic impurity problem is nevertheless qualita- 
tively different for two important reasons. First, as we 
show below, the unstable fixed point separates topologi- 
cally distinct sectors in the Kitaev model, with the strong 
coupling sector associated with non-abelian anyons. Sec- 
ond, the gapless spinons in the Kitaev spin liquid me- 
diate a non-dipolar RKKY interaction proportional to 
SfaSj/s/Rij between distant magnetic impurities pro- 
vided that (a) each impurity couples to more than one 
lattice site on the host and (b) the impurity spin S 7^ 1/2. 
The absence of long-range interaction for 5 = 1/2 impu- 
rities opens a way for local manipulation of the Kitaev 



system. A comparison of Kondo effect and RKKY inter- 
action in graphene [l3, 14 1, a bosonic spin bath [l3| and 
the Kitaev model are shown in Table fl] 

The 5=1/2 Kitaev model d is a honeycomb lat- 
tice of spins with direction-dependent nearest neighbor 
exchange interactions, 



H = -J X J2 °] 
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where the three bonds at each site (see FigO} are la- 
beled as x, y and z. As was shown by Kitaev, the 
flux operators W p = <7f erf 03 04 of erf defined for each 
elementary plaquette p are conserved (see FigQ]), with 
eigenvalues ±1, and form a set of commuting observ- 
ables. Each of the Kitaev spins is represented in terms 
of Majorana fermions bf , tif, b\ , Cj as erf = ibfci, which 
span a larger Fock space, and we restrict to the physi- 
cal Hilbert space of the spins by choosing the gauge [j| 
Di = ibfbfb^Ci = 1. On each a— type bond, ufj = ibfb" 
is also conserved and the ground state manifold corre- 
sponds to a vortex free state where all Wi are equal. 
In the vortex free state, we can fix all Uij = 1 (corre- 
sponds to W p = 1) and the Hamiltonian can be written 
as a theory of non interacting Majorana fermions. The 
reduced Hamiltonian for this ground state manifold is 
given by H = \Y. jk A jk CjC k , where A jk = U a . k if 
j, k are neighboring sites on an a— bond and zero oth- 
erwise. The excited state manifolds (with finite vortic- 
ity) are separated from the ground state manifolds by a 
gap of order J a . Defining the Bravais lattice with a two 
point basis (FigflJ, the Hamiltonian can be diagonalized 
in momentum space, H = | Y, a >o a £ a(q)a t 



2(J X 



^q>0,Q 
,iaq-ni 



q ,a«q,Q> With 



+ J z ) 



e a (q) = ±|/(q)|, /(q) 

and the eigenbasis, a q = c^a + c q ge and a q .i 

Cq.A — c qi .Be~ M ( q ). Here A/B is the site label for the two 
types of sites in Kitaev model, a is the lattice constant 
and a(q) is the phase of /(q). The sum over momenta 
is only over half of first Brillouin zone. e(q) has lin- 
ear dispersion around the Fermi point kp (FiglTJ . For 
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Graphene 


Z2 bosonic spin bath with pseudogap 
density of states p[e) = e . 


Kitaev, honeycomb lattice 


Kondo scaling 


Unstable intermediate coupling 
fixed pt. only for AFM coupling. 
Only AFM flows to strong 
coupling above unstable fixed pt. 


Flow direction is independent of the 
sign of magnetic impurity coupling. 
Unstable intermediate coupling fixed 
pt. for both FM and AFM. 


Scaling same as Z2 bosonic spin 
bath case. However a topological 
transition is associated with the 
unstable fixed point. 


RKKY 


SiaSj a / Rij 







Table I: Comparison of Kondo effect and RKKY interaction in graphene, a Z2 bosonic spin bath with a pseudogap density of 
states and the Kitaev model on the honeycomb lattice. 




Figure 1: (a) Schematic of the Kitaev lattice showing the 
A and B sites and the x, y and z types of bonds, (b) Fig- 
ure showing the reciprocal lattice vectors for the A sublat- 
tice. The Dirac point for the massless Majorana fermions 
is denoted by &f and momentum summations are over the 
(shaded) half Brillouin zone. 

simplicity, and without loss of generality, we henceforth 
assume J x = J y = J z = J. 

Topological Kondo effect - Consider a spin S mag- 
netic impurity coupled to a Kitaev spin at an A site 
(r = 0), 

V K =iY / K a S a b a c A (0), (2) 

a 

for which we perform a standard poor man's scal- 
ing analysis [15j for the Kondo coupling K. Con- 
sider the Lippmann-Schwinger expansion for the 
T— matrix element, {b^\K^ S^b^ c a ,A\{^ a)) (scattering 
of a c— Majorana with momentum q and sublattice in- 
dex a to a &-Majorana), T = + H , in in- 
creasing powers of K. The first correction to the bare 
T— matrix comes from two third order terms (see Figf2|): 
T (3) „ _^ K P S P£i^J2p(KP) 2 (SP) 2 . Here p(e) = 
(l/2irvp)\e\ = C|e| is the density of states and D is the 
band edge energy. If either the impurity is a S = \ spin, 
or the Kondo interaction is rotationally symmetric, the 
above contribution renormalizes the Kondo coupling con- 
stant. However for S 7^ | with anisotropic coupling, new 
terms are generated and one needs to go to higher or- 
der diagrams to obtain the scaling of these new coupling 
terms. 

Just as for the Kondo effect in graphene 13], we also 
need to consider the change in the density of states 
with bandwidth. This gives a contribution K — > 
K(D'/D) r , (£>' = D - \SD\). For S = 1/2 or for sym- 




Figure 2: (a and b) Third order contributions to the 
T— matrix. Site-diagonal scattering corresponds to i = j 
and site off-diagonal scattering, where relevant, corresponds 
to i / j. Thin solid lines correspond to c— Majoranas while 
dashed ones to b— Majoranas. Thick solid lines represent 
the impurity spin, (c and d) New vertices generated by off- 
diagonal scattering. 

metric impurity coupling we thus have 

5K = -K 5 -^- (2K 2 a 2 CDS(S + 1)/J - l) . (3) 

The effective coupling K has an unstable fixed point at 
K c = ^/J/[2a 2 p(D)S(S + 1)] - J/S. Here we used D < 
J and C ~ l/(Ja) 2 . Clearly for K > K c , the coupling 
flows to infinity independent of the nature of coupling 
(ferromagnetic or antiferromagnetic), while for K < K c , 
the coupling flows to zero. 

If the impurity couples to more than one Kitaev spin in 
a plaquette, new contributions arise from site off-diagonal 
scattering (i 7^ j in Fig[2] a, b) . Adding all these contri- 
butions, we find this also leads to a similar Kondo effect 
as was discussed above for the single site coupling case. 
There are also new terms of second order in K that are 
generated (see FigJ5]). The term corresponding to Fig. 

Efc is ~ K^K^S^Sh^. When projected to the vor- 
tex free ground state, it becomes ~ {K^) 2 {S^) 2 gener- 
ating anisotropic potential for the impurity spin. The 
second term, FigHH, is — (K 13 ) 2 {S^) 2 CiCj/J, which, as 
we shall see below, contributes to the long range interac- 
tion among impurity spins. 

A remarkable property of the Kondo effect in Kitaev 
model is that the unstable fixed point is associated with 
a topological transition from the zero flux state to a fi- 
nite flux state. The strong coupling limit amounts to 
studying the Kitaev model with a missing site or cutting 



3 



the three bonds linking this site to the neighbors. Kitaev 
has shown [B| that such states with an odd number of 
cuts are associated with a finite flux, and also that these 
vortices are associated with unpaired Majorana fermions 
and have non-abelian statistics under exchange. Below 
we present another, perhaps more physical, way of ap- 
preciating this result. 

For the Hamiltonian H = Hq + Vk, the three plaque- 
ttes Wi , W2 and W3 that touch the impurity site are no 
longer associated with conserved flux operators, while the 
flux operators that do not include the origin remain con- 
served. The three plaquette operator Wq = W1W2W3 
is still conserved and Wq — 1 in the ground state of 
the unperturbed Kitaev model. The composite opera- 
tors t x = W 2 W 3 S X , T y = W$W\S y and t z = W{W 2 S Z , 
where the S a are Pauli spin matrices corresponding to 
the impurity, are also conserved. The t q 's do not mu- 
tually commute and instead obey the SU(2) algebra, 
[T a ,T P ] = 2ie af)l T"i . This SU{2) symmetry, which is 
exact for all couplings is realized in the spin- 1/2 repre- 
sentation ^(t q ) 2 = ij. Thus all eigenstates, including 
the ground state are doubly degenerate. 

Consider a strong (antiferromagnetic) coupling limit 
Jk 00. The low energy states will be the ones in which 
the spin at the origin forms a singlet |0) with the impu- 
rity spin, = \ipK~) ® |0). \ipK~) now stands for the 
low energy states of the Kitaev model with the spin at 
the origin removed. To sec the action of the SU{2) sym- 
metry generators on these states, we note that they can 
be written as t q = W a <g> a a <g> S a and W a do not in- 
volve the components of the spin at the origin, a a . We 
then have T a \ip) = -(W a \ipK-)) (g> |0). Thus, in the 
strong coupling limit, the symmetry generators act non- 
trivially only in the Kitaev model sector. This implies 
that the low energy states of the Kitaev model with one 
spin removed are all doubly degenerate, with the double 
degeneracy emerging from the Kitaev sector. This im- 
plies there is a zero-energy mode in the single particle 
spectrum. The two degenerate states correspond to the 
zero mode being occupied or unoccupied. The same ar- 
guments for the double degeneracy in the Kitaev sector 
may be repeated for the ferromagnetic strong coupling 
case. 

Let us examine the structure of the zero mode. Remov- 
ing a Kitaev spin creates three unpaired b— Majoranas at 
the neighboring sites, say, 63, b% and b\ . Note that ib%b\ 
is conserved and commutes with all the conserved flux 
operators Wi but not with the two other combinations 
ib\b\ and ib^b^ - then one can choose a gauge where the 
expectation value of ib%b\ is equal to +1 so that these 
two b modes drop out of the physics and we equivalently 
have one unpaired b— Majorana. The unpaired 6| Majo- 
rana has dimension y2; so there must therefore exist an 
unpaired Majorana mode in the c sector (again of dimen- 
sion y/2) so that together these two give the full (doubly 
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Figure 3: Schematic of a typical long-distance impurity in- 
teraction mediated by a pair of propagating c— Majoranas 
(shown by dashed lines) emanating from the ends of a Kitaev 
bond. Another contribution comes from crossed trajectories 
(i,A«/,B' and j, B i',A'). 



degenerate) zero energy mode. We note that while the 
&| mode is sharply localized, the wave function of the c 
mode can be spread out in the lattice. 

That the ground state energy corresponding to a finite 
flux state (pinned to the defect) is lower than the zero 
flux case in the Kitaev model with a missing site has also 
been recently shown numerically fislj . 

RKKY Interactions - In the absence of impurities, 
in the ground state manifold (vortex free state) , we have 
only nearest neighbor Kitaev spin correlations. This is 
because each Kitaev spin is a bilinear of a massless c 
Majorana and a massive b Majorana, and the b Majo- 
ranas have only short range correlations. Suppose we 
now add impurities which may each be locally coupled 
to more than one Kitaev spin. Distant impurities can 
interact only if they are coupled via the massless Ma- 
joranas. By contracting b Majoranas locally by second 
order perturbation in the Kondo coupling, we gener- 
ate vertices of the type (K a S a ) 2 CiCj / J, where i and j 
are not farther than nearest neighbor. Note that since 
cf = 1, these vertices will contain massless Majoranas 
only when i and j belong to different sites. This effec- 
tively means that two distant impurities coupled to a 
single Kitaev site each cannot interact. However when 
the impurities interact with more than one neighboring 
Kitaev spin, we shall see that a long range interaction 
of the spins is possible. As an example, we analyze 
the interaction when the two impurities are at the cen- 
ters of distant hexagons. The interaction term is Vk = 
^ J e^i, a K a S^c ] +iY, ]&c ^ a K a S^c J . The ef- 
fective interaction generated involving c— type Majorana 
fermions at two neighboring sites (i € A, j 6 B) is given 
by (see Fig[2b,d) V cS = (J^) 2 ^) 2 ^. 
Here cty refers to the z— component when neighbor- 
ing sites i and j are along a z— bond, etc. Now 
the interaction between the two impurity spins is 
given by the second order term in V c s (or, equiva- 
lently, fourth order in K). These terms are of the 

type _ £ <(X«« f{ST] ) 2 (K^nS^' f Cl c jCl .c r ) . Per- 
forming the fermionic averaging, the contribution (see 
Fig. to the long range interaction from the pair of 
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bonds ij and i'j' is 



J 12 



4vfJ 2 t^ 3 



1 + cos(2d(k F )) - 2 cos(2k F • R 12 ) 

n X2 



(4) 



Note that for spin-1/2 impurities, (S a ) 2 — 1/4 and for 

isotropic coupling where £bo„d pa irs( S T 3 ) 2 ( S 2*' 3 ' ? = 
const., no long-ranged interaction is generated. 

To summarize, we studied the effect of impurity quan- 
tum spins coupled to the ground state manifold of the 
Kitaev model in the gapless spin-liquid state. We found 
an unusual Kondo effect with an unstable fixed point 
demarcating a topological transition between zero flux 
and finite flux sectors. Where more than one impurity is 
present, we showed that under certain circumstances, the 
massless spinons in the Kitaev model mediate a higher 
order (non-dipolar) RKKY interaction between distant 
impurity spins. The topological transition and the non- 
dipolar impurity interaction make the Kitaev Kondo ef- 
fect qualitatively different from the Kondo effect in some 
bosonic spin liquids that also have an unstable fixed 
point. We expect a similar scaling for the Kondo coupling 
and RKKY interaction for other spin models [17jwith a 
similar Majorana structure and phase diagram. 

In the strong Kondo coupling limit we showed that a 
non-abelian anyon is created consisting of an unpaired 
b— Majorana localized in its vicinity and its delocalized 
c— counterpart. The localized b— Majorana at the finite 
flux defect is very reminiscent of the localized Majoranas 
in the cores of half vortices in p— wave superconductors 
[HI- One difference, as pointed out by Kitaev, is that 
in p— wave superconductors, the currents associated with 
the fluxes are charged (and thus interact with impu- 
rity potentials). Besides, the individual charges acquire 
abelian phases of their own. Another difference is that 
in the Kitaev model, the full (doubly degenerate) zero 
mode is made of a b and a c Majorana while in the su- 
perconductor, two vortex core Majoranas make a zero 
mode 23]. p— wave paired ground states in the Kitaev 
model and its vortex excitations have also been studied 
inRef. [H. 

We note that if we could adiabatically move the im- 
purity spin (say by making it on a STM tip), then we 
could move the anyons and thus perform the braiding 
operations as is required for quantum computation. This 
method is more practical than perturbing an entire pla- 
quette with terms like fjWi that involve at least six Ki- 
taev spins. 

The temperature dependent impurity susceptibility 
can be measured by NMR Knight shift experiments 
[20}. For weak Kondo coupling, the magnetic suscep- 
tibility of the impurity is Curie-like with logarithmic cor- 
rections, which flow to zero regardless of the sign of 
Kondo coupling. In the strong coupling limit, we con- 
sider antiferromagnetic coupling for S — 1/2 impuri- 



ties, where the anyons are vacancies in the Kitaev lat- 
tice. The magnetic susceptibility has been shown 16 1 
to have the form x(T) ~ T ln ^ D yj* for a pair of nearby 
vacancies on the same sublattice while for a single vor- 
tex, x{T) ~ \n(D/T). In the absence of vacancies, the 
low temperature magnetic susceptibility is small because 
of the spin gap in the Kitaev ground state. The nu- 
clear relaxation rate T^ 1 may also be used as an im- 
purity probe [2lj]. Using x(T,lu) ps (x _1 (T,0) - itu)- 1 
and Tf 1 cx A 2 I(I + l)T[Imx(T,w)/w] w _ > Q ) where A is 
the hyperfine coupling of spin—/ nuclei with the defect, 
one gets T^ 1 ~ T[ln(D /T)] 2 for the single vortex and 
Tf 1 - l/T[ln(£)/r)] 2 for the two nearby defects case, 
both of which qualitatively differ from the Korringa law 
T^ 1 r-j T for the usual Kondo effect at low temperatures. 
There are already numerous proposals in the literature 
how a Kitaev model could be realized [22], and we are 
hopeful that eventually these novel impurity effects may 
also be experimentally studied. 
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